Fourier processing of quantum light 
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It is shown that a classical optical Fourier processor can be used for the shaping of quantum 
correlations between two or more photons, and the class of Fourier masks applicable in the multi- 
photon Fourier space is identified. This concept is experimentally demonstrated using two types of 
periodic phase masks illuminated with path-entangled photon pairs, a highly non-classical state of 
light. Applied first were sinusoidal phase masks, emulating two-particle quantum walk on a periodic 
lattice, yielding intricate correlation patterns with various spatial bunching and anti-bunching effects 
depending on the initial state. Then, a periodic Zernike-like filter was applied on top of the sinusoidal 
phase masks. Using this filter, phase information lost in the original correlation measurements was 
retrieved. 

PACS numbers: 42.30.Kq, 42.50.Dv, 03.67.Ac, 03.67.Bg 



Fourier processing is a well-established method in sig- 
nal processing, where a transformation is applied to the 
signal in the reciprocal Fourier domain. In classical op- 
tics, a lens converts an optical field in the back focal plane 
to its spatial Fourier transform at the front focal plane, 
where optical operations can be applied to directly ma- 
nipulate the signal in the fourier domain This basic 
scheme has been used for the control and manipulation of 
the light intensity distribution in a number of optical sig- 
nal processing methods, from matched filtering to phase- 
contrast microscopy 1]. One of the most frequently used 
setups is that of the 4-f filter, where a spatial mask at the 
common focal plane of two lenses performs the Fourier 
processing of the input signal (see inset in Fig. [T]). 

In recent years, there is a growing interest in the gen- 
eration, control and detection of non-classical, quantum 
light [2[. This is mostly due to its exciting applications 
in quantum information processing [3[ , quantum metrol- 
ogy 13, quantum imaging [Bj], and quantum lithogra- 
phy [6(. In contrast to classical light, many properties of 
quantum light cannot be revealed through first-order cor- 
relations, such as the light intensity, and require the mea- 
surement of higher order correlations 0, 0] ■ Since optical 
Fourier processing proved to be a powerful tool for con- 
trolling the intensity distribution of classical light, exam- 
ining the possibility of using this technique for the con- 
trol of quantum correlations between photons seems very 
promising. This approach is further motivated by recent 
works applying Fourier optics to two photon states [§- 
13j| . where it was shown that the optical Fourier trans- 
form of two photons in one spatial dimension is given by 
the two-dimensional Fourier transform of the two-photon 
wavefunction [id lll| . 

In this work, we first establish an analogy between the 
Fourier processing of classical light (intensity) to that 
of quantum light (correlations) by studying the effect of 
the Fourier-plane mask on photon correlations. We show 
that the effect of a one-dimensional mask in the space of 
the two-photon wavefunction is equivalent to the two- 



dimensional, external product of the one-dimensional 
mask with itself. This defines the class of Fourier fil- 
ters that can be applied in the two-photon correlation 
space using a classical, one dimensional mask. We then 
use these results to experimentally demonstrate the shap- 
ing of non-classical correlations between two photons by 
spatial Fourier processing with various masks. 

For the experimental demonstration, path-entangled 
photon pairs, a highly non-classical state of light, were 
introduced into a 4-f filter. The resulting two-photon cor- 
relation maps were measured for two different types of 
phase masks implemented at the Fourier plane by a spa- 
tial light modulator (SLM) . The first mask extended the 
correlation object by effectively implementing quantum 
walk of two particles and the second mask was 

used as a Zernike-like filter, similar to that used in classi- 
cal phase-contrast microscopy, for retrieving two-photon 
phase information. 

Before we discuss the details of the experiments and 
the experimental results, we would first like to review 
the application of the optical Fourier transform to more 
than one photon, and use it to establish the analogy be- 
tween single-photon and multi-photon Fourier process- 
ing. For the simplicity of the following discussion we 
consider monochromatic light, omit the polarization de- 
gree of freedom, and consider only polarization main- 
taining, linear optical devices. A single-photon tran- 
sition operator describing such a device has the form, 
U=J dx dxiU(x 0l Xi)w (x )a(xi). The operator d{xi) 
\w(x )] annihilates [creates] a photon at the coordinate 
Xi [x ] of the input [output] plane of the device, and 
U(x ,Xi) is the (complex) amplitude of that transition. 
The coordinates can be either one- or two-dimensional, 
depending on the studied system, and the operators obey 
bosonic commutation relations. We then introduce at the 
input plane of the optical device a two-photon state, 

*m = J dx l idxi 2 As{x t i,Xi2)a){x l 2)a\xii)\Q), (1) 

where Asixn , Xi2)—As{xi2, xn) is an exchange- 
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symmetric complex probability amplitude. The spatial 
intensity-correlation function of this state is given by 
r(xji, Xi2) = \Ag{xn, Xi2)\ 2 ■ For two non-interacting 
photons passing together through the device, the output 
two-photon state is obtained from the input state by 
applying the single-photon transition operator twice, 
^out—U 2 ^>in- One then obtains, 

$out = / dx ol dx o2 Bs(x ol ,x O 2)b' !! (x O 2)P(x ol )\0), (2) 



where the output complex probability amplitude, 
B s (x i,x o2 )=B s (x o2 ,x i), is given by, 



B s (xoi,x o2 ) = J dXiidx i2 As{x i i,x i 2) ■ (3) 
■U(x o2 ,x i2 )U(x ol ,xa), 

The two-photon transition amplitude is thus just the 
external product, U(x n2 , Xg)U(x n i , Xn), of the single- 
photon transition amplitude with itself [l0( • In the same 
manner it can be shown that this property holds for an 
arbitrary number of non-interacting photons. 

We note that this result, though derived here for pho- 
tons, which are bosons, applies also to fermions, provid- 
ing that an exchange- anii-symmetric initial probability 
amplitude is being used. 

We now apply this general result to the components 
of the Fourier processor. We begin by treating the 
lenses, which perform the Fourier transforms from the 
input to the Fourier plane, and from there to the out- 
put. Eqs. ([2]) -([3]) applied to a one dimensional Fourier 
transform, like that performed, for example, by the first 
lens, ULi(xf,Xi)—e~ 27rlXfXi ^ x ^ 1 , recover the known result 
that in the two photon coordinate space this transforma- 
tion is equivalent to a two-dimensional Fourier transform, 
^i(x fu x n )U L1 {xf 2 ,x i2 )=e- 2 ™^ x "+ x f* x ^/ x h 0, 
1 1 11 ] . Here xj is the coordinate in the Fourier plane, /i 
is the focal length of the lens, and A is the wavelength 
of the incident light. This result applies to the second 
lens as well. It can also be generalized to any number 
of photons, N, and spatial dimensions, D, since the ex- 
ternal product of N D-dimensional Fourier transforms is 
just the Nx D-dimensional Fourier transform. 

The second component in Fourier processing is the 
mask implemented at the Fourier plane of the 4-f fil- 
ter. Since the action of a mask with complex ampli- 
tude transmission M{xf) on a single-photon probabil- 
ity amplitude is A{xj) -¥ M(xf)A(xf), from Eq. ((3|) it 
follows that its action on a two-photon probability am- 
plitude is A s (xfi,Xf 2 ) -> M(xfi)M(xf 2 )A s (xfi,Xf 2 ). 
This represents a mask in the two-dimensional Fourier 
plane, {2/1,2/2}, which is the external product of the 
one-dimensional mask M(x/) with itself. This result re- 
flects the fact that a single linear system can act only 
locally on the photons, and induces the same transfor- 
mation to each photon. Since unlike the two-dimensional 



Fourier transform, a general two-dimensional mask is not 
necessarily a product of a one-dimensional transforma- 
tion by itself, only this special class of two-photon masks 
can be implemented using a classical, linear mask. The 
mask, therefore, is the element that limits the ways in 
which a classical, linear Fourier processor can process a 
non-classical two-photon wave-function. The implemen- 
tation of a general mask would require the two photons to 
interact on the mask. This is not yet feasible with cur- 
rent technology, though progress is being made in this 
direction [17] . Nevertheless, even with linear masks, sev- 
eral important processing algorithms can still be imple- 
mented. For example, a Fourier mask transmitting only 
in a small aperture will block fermion pairs, while boson 
pairs will go through. This is since such a mask is equiv- 
alent to a small square aperture on the main diagonal of 
the two-particle Fourier space, and thus transmits only 
the zero order with respect to particle exchange. In con- 
trast to boson pairs, fermion pairs, which wavefunction is 
cxchange-anti-symmetric, do not have such a zero order. 

For the experimental demonstration of two-photon 
Fourier processing we use phase-only masks. In contrast 
to the aperture mask discussed above, phase-only masks 
have, in principle, no loss. On the other hand, they can- 
not change the total two-photon flux. They can, how- 
ever, controllably redistribute it among different output 
two-photon configurations, as we demonstrate below. 




FIG. 1: Experimental setup. Polarized light from a laser 
diode (LD) is split into two beams by a half-wave plate (HWP) 
and a polarizing beam splitter (PBS). The distance between 
the beams is controlled by a movable retroreflector (RR), and 
their relative phase is controlled by a movable mirror (M) on 
a piezoelectric transducer (PZT). The beams undergo down- 
conversion in a nonlinear crystal (BBO). The pump beams 
are diverted by a dichroic mirror (DM) and interfered on a 
camera (CCD). The interference pattern is used together with 
the PZT for actively stabilizing the relative phase between the 
beams. The down-converted light goes through a 4-f filter 
equipped with a spatial light modulator (SLM), and finally 
detected by two fiber-coupled detectors (Dl, D2). Coinci- 
dence detection is manifested by a coincidence counting unit 
(CCU). SF is a spatial filter, QWP is a quarter- wave plate, 
and BS is a beam-splitter. Purple lines represent the pump 
beams, and red lines represent the down-converted light. The 
inset shows a scheme of a 4-f filter. 

The experimental setup is described in Fig. [1] 
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We first generate a two-path-entangled, two-photon 
state of the form, ^(|2 Q , b )+e^|0 a , 2 b )), that is, 

As(xii,x i2 )=^[p{x l i-x a )p(x t 2-x a )+e t ' f 'p(x il -x b )p(x t 2- 
x b )], where x a (x b ) is the central lateral position of 
the first (second) path, and p(xi) is its lateral pro- 
file. This is achieved by down-conversion of a split, 
phase-stabilized laser beam [l8j]. The down converted 
light then goes through a 4-f filter equipped with a 
one-dimensional, phase-only SLM at its Fourier plane. 
Using the SLM, various one-dimensional phase masks 
can be created. In this work we use only periodic phase 
patterns, and set their period to match the separation 
between the two initial modes, such that the output 
intensity pattern is also in the form of discrete modes 
[see Fig. EJa)]. Finally, the total output intensity 
distribution, as well as the output intensity correlations 
between any two modes q and r, T q r , are recorded 
by two fiber-coupled single-photon detectors. Here 
r 9,<-= j™ w f™ w dx idx o2 \Bs{x i-x q ,x 2-x r )\ 2 7 where the 
mode q (r) is centered about x q (x r ), and 2w is the 
width of the collecting fibers. 
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FIG. 2: (a) Schematic of the beam shaping system with a sinu- 
soidal phase mask applied on the SLM, showing beam prop- 
agation and resulting intensity patterns, (b)-(c) Measured 
output intensity patterns for different values of the phase am- 
plitude, A p , for, (b), a single input beam, and (c), two beams 
of path-entangled photon pairs. 

The intensity distribution of the initial state is spa- 
tially limited to two spots. This is true also for the cor- 
responding correlation pattern, since due to the bunched 
nature of the photons in the initial state, the correla- 
tion pattern is composed of two spots on the main diag- 
onal. In order to demonstrate correlation manipulation 
similar to Fourier image processing, we first use a 4-f fil- 
ter with a sinusoidal phase mask, M ^ X f)—e lApCOs(27TUXf \ 
where A p is the phase mask amplitude, and v is the spa- 
tial frequency of the modulation. This filter implements 
the tight-binding model 0, Qjl] , and emulates the quan- 
tum walk [l4l | of the two entangled photons on a pe- 



riodic lattice. This is since the phase mask is applied 
in the Fourier space - the space of the eigen-modes of 
any periodic system, and since the induced phase shift 
is proportional to the sinusoidal tight-binding dispersion 
relation [20j. This filter thus results in spread-out in- 
tensity distributions [l, 14, 19[, and intricate correlation 
patterns that dep end on the phase, 0, between the two 
initial modes [15| . 

Indeed, the measured output intensity patterns for 
a single input beam [two input beams], shown in 
Fig. [2Ib)[(c)], are practically identical to those observed 
in other optical systems implementing the tight-binding 
model, such as a periodic lattice of evanescently coupled 
waveguides 
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The incident light, originally con- 
centrated in one or two input modes, spreads among a 
number of output modes increasing linearly with A p . 




FIG. 3: Experimental observation of non-classical correla- 
tions using a sinusoidal phase mask, (a) Second-order spatial 
correlation function, r q , r , of the input state, as measured at 
the output when no mask was applied. The black frames 
represent the expected values, (b) A section of the phase pat- 
tern created by the applied sinusoidal phase mask in the two- 
photon Fourier plane, (c)-(f) Resulting correlation patterns 
for different phases, (j>, of the input states. Anti-bunching is 
observed for 4> — (c), while for <j> = ir (d) bunching is ob- 
served. For the two intermediate cases, (j> — —tt/2 (e), and 
(/> — tt/2 (f), both effects are seen simultaneously. 
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For the correlation measurements we apply a sinusoidal 
phase mask with A p = 0.867T. The corresponding inten- 
sity pattern is indicated by the dashed line in Fig. [2fc) . 
The phase of the two-dimensional mask that results in the 
two-photon Fourier space is presented in Fig. (3][b) . The 
resulting correlation maps are presented in Figs. EHc)- 
[3jf), for four different initial phases, </>. The full bars 
present the measurement, while the empty bars present 
the theoretical prediction. For comparison, the correla- 
tion map of the initial state, as measured at the output 
when no mask was applied, is presented in Fig. G^a). 

It is seen that, depending on 0, the two photons, which 
are initially bunched together, can emerge anti-bunched 
[4>=0, Fig. OUc)], bunched [</>=7r, Fig. E^d)], or in any 
super-position of these two extremes. Of these super- 
positions, those generated for the two cases of (f>=— tt/2 
and <^—tt/2 have identical correlation maps [See Figs.[3{e) 
and[3^f)]. However, since the initial phase is different 
between the two cases, so should be the phases in the final 
patterns, and we would therefore like to have a procedure 
that will reveal this hidden information. 

The second Fourier processing algorithm, which will be 
performed on the obtained correlation patterns, would 
thus be dedicated to this purpose. We therefore add on 
top of the sinusoidal phase grating a second phase grat- 
ing, implementing an algorithm designated to reveal the 
phase information lost in the correlation measurements 
discussed above. We note that this is equivalent of build- 
ing a second 4-f filter and adding the second grating at 
its Fourier plane, since then the two Fourier planes would 
just be imaged one on top of the other. 

The specific grating we use for this purpose is inspired 
by phase-contrast microscopy. In this method, the zero- 
order of the image is phase-shifted by tt/2 with respect to 
all other orders using a designated Fourier phase filter, 
known as the Zernike filter. For a transparent, appar- 
ently uniform object, with a small spatial phase-shift dis- 
tribution, this leads to interference of the non-zero orders 
with the zero-order, and the creation of intensity fluctu- 
ations proportional to the local phase-shifts For a 
periodic system, the zero-order is at the center of each 
reciprocal lattice unit cell. In order to approximate the 
Zernike filter for the case of two photons, we applied a 
7r/4 phase to the central quarter of each phase lattice unit 
cell (with respect to the original sinusoidal phase). In the 
two-dimensional space of the two photon wavefunction 
this filter shifts the central part of each two-dimensional 
unit cell by tt/2. There is, however, a byproduct: a 7r/4 
shift along the main axes of each cell. The filter phase 
is shown in Fig. |3£a). Our filter thus only approximates 
the Zernike filter. Furthermore, the initial correlation 
patters [Figs. E^e) and [3jf )] are non-uniform, and the 
predicted phase shifts are not small. Nevertheless, the 
filter can still recover qualitative phase information, and 
differentiate between the two cases. 

This is shown in Fig. @] First note that, as seen in 
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FIG. 4: Retrival of phase information using a Zernike-like 
filter, (a) A section of the additional phase pattern applied 
on top of the sinusoidal phase [Fig. EJb)]. (b) A compari- 
son between the output intensity distribution with (dark red) 
and without (light blue) the additional phase. The number 
of each peak is indicated next to it. (c)-(f) Resulting corre- 
lation patterns for the same initial phases as in Fig. EJc)-(f). 
While for 0=0 (c) and <f>=Tv (d) no qualitative change is seen, 
the two intermediate states, 4>=— tt/2 (e) and 0=7r/2 (f) now 
show different and opposite correlation patterns, revealing the 
phase difference between them. 

Fig. QJb), the filter has almost no influence on the in- 
tensity distribution. This is due to its unitarity and its 
periodicity, matching that of the sinusoidal phase mask. 
The filter also does not qualitatively change the correla- 
tion patterns for the cases of (f>=0 [Fig. 0|c)] and tp—ir 
[Fig. QJd)] . This is since the wavefunctions in these cases 
are all real, and thus the Zernike-like filter has little or 
no influence on them. In sharp contrast, for the cases 
of 0=-tt/2 [Fig. gfe)] and 0=tt/2 [Fig. gff)], the ap- 
plication of the Zernike-like filter changes the correlation 
patterns, and creates a clear difference between these two 
cases, which now show opposite, asymmetric correlation 
patterns. 

In conclusion, we have demonstrated that the correla- 
tions between path-entangled, two photon states can be 
manipulated using a classical Fourier processor. Though 
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the applicable masks are limited to local operations, 
many types of manipulations are still possible. We used a 
phase-only SLM, and applied phase masks both to create 
an elaborate correlation pattern through the emulation of 
quantum walk on a periodic lattice, and to recover hid- 
den phase information. For doing so we used methods 
adapted from classical Fourier processing, such as phase 
contrast microscopy. The analogy between image and 
correlation Fourier processing can be further exploited by 
using filters that change the amplitude of the light as well 
as the phase. Examples of such filters include low-, high-, 
or band-pass filters, the schlieren filter, and matched fil- 
ters. Furthermore, additional degrees of freedom can be 
utilized by considering two-dimensional and/or polariza- 
tion changing masks. Moreover, this analogy is applica- 
ble also in the time-domain, where it can be applied for 
the tailoring of temporal correlations between photons 
through pulse-shaping techniques [21|. We thus believe 



that the use of Fourier processing methods may open new 
avenues for the manipulation of non-classical light. 
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